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lepton number, focusing on the L-conserving contribution to the flavoured CP asymmetries. 
We find that, contrary to previous claims, successful leptogenesis is feasible for masses of 
the lightest heavy neutrino as low as M\ ~ 10 6 GeV, without resorting to the resonant 
enhancement of the CP asymmetry for strongly degenerate heavy neutrinos. This lower 
limit renders thermal leptogenesis compatible with the gravitino bound in supersymmetric 
scenarios. 
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1. Introduction 

Two of the main evidences of physics beyond the Standard Model (SM), namely the ob- 
served (tiny) neutrino masses and the baryon asymmetry of the universe (BAU), are natu- 
rally explained in the context of the seesaw mechanism [1, 2, 3, 4]. In type I seesaw models, 
the SM is extended with at least two singlet Majorana neutrinos which can produce the 
observed BAU via leptogenesis [5, 6]: a lepton asymmetry is dynamically generated in the 
out of equilibrium decay of the heavy Majorana neutrinos, and then partially converted into 
a baryon asymmetry due to (B + L)-violating non-perturbative sphaleron interactions [7]. 

For a very hierarchical spectrum of heavy singlet neutrinos M\ <C M2 <C M3, the 
L-violating CP asymmetry generated in the decay of the lightest singlet has an upper 
bound proportional to Mi, the so-called Davidson-Ibarra (DI) bound [8]. This implies a 
lower bound ~ 10 9 GeV for the mass of the sterile neutrinos in order for iVi-dominated 
leptogenesis to be successful. Careful numerical studies show that the DI bound can be 
evaded for moderate hierarchies, e.g. the lower bound on Mi is relaxed by more than one 
order of magnitude with respect to the hierarchical limit one for M3/M2 ~ M2/M1 ~ 10 [9]. 
However to reach these low values of Mi some unlikely cancellations are needed, which are 
not motivated by any underlying symmetry. Flavour effects [10, 11, 12, 13, 14, 15] do not 
substantially change this result. The authors of [16] make an analysis of the parameter 
space for successful leptogenesis relaxing also the condition of hierarchical neutrinos and 
including the L-conserving part of the flavoured CP asymmetries, finding a lower bound 
Mi > 10 s GeV. 

Such lower bound on Mi in turn yields a lower bound for the reheating temperature, 
Trh, of the same order, since to thermally produce the neutrinos Mi < 5Tr# [17, 18, 19]. 
On the other hand, in supersymmetric scenarios gravitinos are copiously produced in the 
high temperature plasma, and their late decay can jeopardize successful nucleosynthesis 
(BBN), leading to an upper bound on Trh, which depends on the gravitino mass [20, 21]. 
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If the gravitino is unstable, for gravitino masses, 7713/2, in the natural range from 100 GeV 
to 1 TeV, and within the minimal supergravity framework, Trh should be smaller than 
10 5 - 10 7 GeV, while for m 3/2 > 10 TeV, T RH can be of order 10 9 - 10 10 GeV [22]. As 
a consequence, in supersymmetric thermal leptogenesis there is some conflict between the 
gravitino bound on the reheat temperature and the thermal production of heavy neutrinos. 

There are several possible ways out of this conflict: for instance, if the gravitino is 
stable, the nucleosynthesis bound depends on the next-to-lightest supersymmetric particle, 
but Trh ^ 10 9 GeV can be obtained for m 3 / 2 ^ 10 GeV [22]. Also, one can avoid 
the DI bound resorting to resonant leptogenesis, i.e., a resonant enhancement of the CP 
asymmetry which occurs when there are at least two strongly degenerated heavy neutrinos, 
such that M2 — M1 ~ Ttv, being Fn their decay width [23, 24]. In this scenario, leptogenesis 
is feasible at much lower temperatures, T ~ 0(1 TeV) [25, 26, 27]. 

Although such quasi-degeneracy of the heavy neutrinos might seem unnatural, there 
are well motivated seesaw models which yield a heavy neutrino quasi-degenerate spectrum, 
namely those with an approximately conserved B — L [28, 29]. In these models the tiny 
neutrino masses are proportional to small lepton number-breaking parameters, which are 
technically natural since a larger symmetry is realized when they vanish [30]. Moreover, the 
heavy neutrinos can be much lighter than in the generic seesaw, within the energy reach of 
LHC, and there can be a large active-sterile neutrino mixing. Also, lepton flavour violation 
rare decays as well as non-unitarity of the leptonic mixing matrix are present even in the 
limit of conserved B — L, and therefore they are unsuppressed by the light neutrino masses 
[31, 32, 29]. As a consequence, much attention has been devoted recently to this class of 
low scale seesaw models, since they have a rich phenomenology both at LHC [33, 34, 35] 
and at low energy charged lepton rare decay experiments, such as [i — > cy, and also lead 
to successful resonant leptogenesis [36, 37, 38]. 

In order to see more clearly the structure of seesaw models with small violation of 
B — L, it is useful to expand the low energy effective Lagrangian as 




where the dimensionless couplings c d=5 ,cf =& are assumed to be of 0(1). 

The only operator of d = 5 is Weinberg's operator, responsible for neutrino masses, 
which is suppressed by the lepton number violating scale A.ln, while flavour changing but 
lepton number conserving d = 6 operators Of =6 are suppressed by a different scale Afl- 
In models with approximately conserved lepton number, there is a separation of scales: 
A-fl ~ 0(TeV) can be related to the heavy neutrinos mass scale, but Aln 3> Afl, since 
light neutrino masses only appear when small L-violating perturbations in the Yukawa 
couplings and/or the singlet neutrino mass matrix are introduced; in this scenario, Aln 
does not correspond to any particle mass. 

It has been noticed that even if the heavy neutrino that generates the BAU is not 
quasi-Dirac, or the mass splitting is outside the resonant regime, in seesaw models with 
almost conserved B — L the scale of leptogenesis can be lower than in the standard seesaw 
[39], provided flavour effects are at work. This is so because in these models leptogenesis 
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can be driven by the d = 6 operators, that contribute only to the L-conserving part of 
the flavoured CP-asymmetries which escape the DI bound, even if the heavy neutrinos are 
hierarchical, because these operators are not linked to neutrino masses. The drawback is 
that the same d = 6 operators induce large lepton flavour violating (but lepton number 
conserving) processes, which tend to equilibrate the flavoured asymmetries diminishing the 
total lepton asymmetry [40, 41]. 

In this paper we reanalyze the possibility of having successful leptogenesis driven by the 
purely flavoured L-conserving contribution to the CP asymmetries, in the context of seesaw 
models with small violation of B — L. In the numerical analysis of [16] such contribution 
was included, however models with almost conserved B — L (which can lead to larger CP 
asymmetries) were not considered in detail and the flavour equilibration processes were not 
taken into account. A study of the lower bound on M\ for approximately B — L-conserving 
models was carried out in [39], but their cross sections are different from ours and the 
analysis of the parameter space seems to be non-exhaustive. We have computed the cross 
sections for the crucial lepton flavour changing processes which lead to flavour equilibration, 
and we have found important differences with respect to previous results [26] . We have also 
taken into account decays and inverse decays of the next-to- lightest neutrinos, neglected in 
[39], and thoroughly scanned the parameter space to find the lowest M\ able to generate 
the BAU within this framework. Moreover, we have included a detailed discussion about 
leptogenesis in models with almost conserved B — L, which depends on the heavy neutrino 
spectrum. 

The outline of this paper is as follows. In Section 2 we discuss the main features of 
leptogenesis within the framework of seesaw models with approximately conserved B — L. 
In Section 3 we write the set of Boltzmann equations (BE) relevant for leptogenesis in this 
scenario and compute the cross sections for the lepton flavour violating (but total lepton 
number conserving) processes, which play a crucial role in models with null (or negligible) 
total CP asymmetry. In Section 4 we perform a detailed analysis of the parameter space 
which leads to successful leptogenesis with only L-conserving CP asymmetries, and we 
conclude in Section 5. 

2. Leptogenesis in models with small violation of B — L 

In this section we want to determine what can be distinct, regarding leptogenesis, among 
different models with small violation of B — L. Let us call N\ the SM fermion singlet that 
is mainly responsible for the generation of the lepton asymmetry. Additionally we call N2 
the SM fermion singlet that makes the most important (non resonant) virtual contribution 
to the CP asymmetry in the N\ decays. They have Yukawa interactions with the lepton 
doublets of the SM l a (a = e, fj,, r) described by the Lagrangian 

£ Y = -\ m rtP^N i e a -\* m £ a P R N, i h. (2.1) 

If B — L is only slightly violated, then for each iVj the conditions (i) or (ii) described below 
must be satisfied: 
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(i) Ni is a Majorana neutrino with two degrees of freedom, whose Yukawa interactions 
violate lepton number and therefore the couplings A a j must be small. This is analo- 
gous to the standard seesaw, so the Aj contribution to the d = 5 Weinberg operator 
is given by 

( c m 5 )«/3 = AciA^j . (2.2) 

(ii) The Aj is a Dirac or quasi Dirac neutrino with four degrees of freedom; this means 
that there are two Majorana neutrinos iV^ and Nu with masses Mj + Hi and Mj — Hi 
respectively. The parameter Hi <C Mj measures the amount of B — L violation, so 
that if B — L is conserved, Hi = and Ni = (Aj/j + iNu)/y/2 is a Dirac fermion. The 
Yukawa interactions can be expressed as 



r \ Z\ d N ih + iNu . ~ t N ih -iNu 

c Y m = -Ki h l P R -j= £ a - \ ai W P R -j= £ a + h.c, (2.3) 

where X' ai <C 1. The terms proportional to X' ai induce lepton number violation even 
when \ii — > and hence they are similar in nature to the ones described in (i). 
Neglecting these terms, the Yukawa couplings of Nih and Nu are equal to \ a ih = ^ 
and X a ii = i^A, respectively. These \ a i can be large, because they do not vanish in 
the B — L conserved limit: in the absence of fa an d X' ai , a perturbatively conserved 
lepton number can be defined, by assigning Ln = 1 to Ni, and L^ a = 1 to the SM 
leptons. 

At leading order in the small B—L breaking parameters (A^j, Hi/ Mi), the contribution 
of a quasi Dirac heavy neutrino to the Weinberg operator is 

- (A - - M~i ~M~i ^ + Xm w^ ~ Wi 13 + • • • (2 - 4) 

As a result, one can reproduce the tiny light neutrino masses with large Yukawa 
couplings \ a i and heavy neutrino masses Mj as low as the TeV scale. Moreover, the 
admixture among singlet and doublet SU (2) neutral leptons (and the corresponding 
violation of unitarity in the light neutrino sector) is of order X ai v/Mi and can be 
large. 

In this work we will consider the cases 

- (iia) Hi < ^N ih ,^Nu (Dirac limit), and 

- (iib) rjv jh ,rjv j( <C Hi ^ Mj (Majorana limit), 

where the decay widths T s are given by 

_ Mj + Hj (At A) _ Mj - Hi (A f A)^ _ 
Nlh 8tt 2 ~ 8vr 2 ' { ' ' 

Instead we leave the study of Hi ~ ^N ih u for future research because it is not clear 
that this case can be described with the simple BE used here [42, 43]. 
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Then, the study of leptogenesis for a large class of models with small violation of B — L 
can be covered by considering the different combinations of options for both N\ and A^. 
Next we comment on some key points for the different possibilities, while the complete 
quantitative analysis is made in Sec. 4. We assume that Aq is lighter than N2, leaving 
some remarks of the opposite case M2 < Mi for Sec. 4. 

I. Aq and N2 satisfy (i): This case is not very interesting since the CP asymmetries e a \ 
in the decays of Aq into leptons of flavour a, e a i = ^^tw^j^ttzww^t^Ia 1 being 
proportional to the square of the Yukawa couplings of N2, are very small, so it is not 
different from the standard seesaw. 

II. Aq satisfies (i) and N 2 (h): The L- violating part of the CP asymmetry, ef 1; is sup- 
pressed by fi2 while the L-conserving part, e^ l5 is not. To see this, let us write 
eai = 4i + as [44, 45] 



e al = ^2 /(aj)Im A* J -A Q i(A t A) J -i + ^ g(aj)lm A* J -A Q i(A t A)i J - 

j=2h,2l j=2h,2l 



(2.6) 



where f(a,j) and g(cij) are functions of cij = Mj/Mf and contain the factor l/(A^A)n. 
The terms proportional to / (a?) and g(a,j) come from the L- violating and L-conserving 
contributions, respectively. To lowest order in ^2, /(02/1) = fip-2i) and g(ci2h) = 
g{a,2i)- Taking into account the alignment between the Yukawa couplings of N21 and 
^2hi A a 2i = a ^a2h with a = i, it is clear that 



eai ->/(a 2 / l )Im A* 2/ l A a i(A t A) 2 / l i (l + a* 2 )+ 

M2-^0 L J 



5(a 2/l )Im A* 2h A a i(A t A)i2/ l 



>1 + |a| 



(2.7) 



+ 25(a 2fl )Im A* 2 / l ^ai(A t A)i2/ l 



Hence we see that the contributions from N2h and N21 add up in the L-conserving 
part and cancel in the L-violating one. 

While the total CP asymmetry e\ = ^ Q e a \ is proportional to fi2 because J2 a e ai = 
the amount of matter-antimatter asymmetry produced during leptogenesis generally 
is not suppressed by \ii due to flavour effects (it is possible to have successful lepto- 
genesis even with ei = [13]). Nevertheless, this asymmetry does vanish or becomes 
proportional to fi2 when the L-violating parameter (A'A)n — > 0, independently of 
whether the density of Aq at the onset of leptogenesis is null or thermal. This is 
because for an initial null density of heavy neutrinos and small values of (A^A)n the 
final baryon asymmetry is proportional to (A^A)^ [19], while for an initial thermal 
Aq density the washouts become negligible and so do the flavour effects, in which 
case the total asymmetry generated is proportional to e± oc fi2- 

III. Aq satisfies (ii) and N2 (i): Since the Yukawa couplings of N2 are small, the virtual 
contribution of N2 to the CP asymmetries e a n and € a \h i n the decays of Nu and 
Nih, respectively, is small. More interesting is the virtual contribution of Nu to 
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e a ih and of Nih to e a u, due to the resonant enhancement of the CP asymmetry 

Y 

for degenerate neutrinos. However these contributions are suppressed by — ^== 

because of the alignment of the Yukawa couplings of Nu and Aq^ for A^ = 0, 
therefore the CP asymmetry cannot reach the maximum value 1/2 of more generic 
resonant leptogenesis models. A detailed study of this case was made in [36], in the 
limit \i\ » ^N u ,N lh - they found that it is possible to lower the scale of leptogenesis 
provided that the parameter which controls lepton number violation in the Yukawa 
couplings e = A'/A is much larger than the one describing the mass splitting between 
the quasi-Dirac heavy neutrinos, eu = AM/M. For instance, in order to obtain 
successful leptogenesis with M = 10 6 GeV (1 TeV), one needs e ~ 1CT 3 and eu ~ 
1CT 8 (1CT 11 ). The case \i\ < T^ lu N lh has been considered in [37], however it has been 
shown that in this maximal resonant regime the Boltzmann picture breaks down and 
the Schwinger-Keldysh/Kadanoff-Baym formalism is required [42, 43]. 

IV. Aq and N2 satisfy (ii): Besides the resonant contributions which are the same as in 
III, now also the virtual contribution of N2 to the CP asymmetries e a u and e a ±h can 
be large, so we focus on these terms. Neglecting the X' ai <C \ a i couplings, the CP 
asymmetries in the decays of Nu and N\h become equal. Here it is very important to 
distinguish between the cases fii 3> Tn u ,n 1h and n\ <C ^N u ,N lh - In the first one, Nu 
and Nih behave as two independent Majorana neutrinos regarding the generation 
of the lepton asymmetry, which can hence be roughly double with respect to II. 
However in the second case Ni is (or effectively behaves as) a Dirac neutrino, i.e. 
lepton number is conserved in its decay, and therefore the only possibilities to end 
up with a non-zero baryon asymmetry is to have important washouts from the two 
Majorana components of N2 (if fJ-2 3> ^N 2i ,N 2h ) or let the sphalerons freeze out during 
leptogenesis [46]. 

Summarizing, in models with small violation of B — L the CP asymmetries in the 
decay of N\ can be enhanced in cases II, III and IV. Since the resonant contribution to the 
L-violating CP asymmetry has been widely studied, we develop in Sec. 4 a quantitative 
analysis which covers all the non-resonant interesting cases, i.e., II and IV with fj,i S> 
^N u ,N lh - With respect to N2, we consider two possibilities, fj,2 3> TN 2l ,N 2h and [12 
TN 2t ,N 2h (see [46] for leptogenesis with both N\ and N2 satisfying fj,i <^ ^N u ,N ih )- 

3. Boltzmann equations 

Motivated by the above discussion we consider a scenario for leptogenesis involving three 
fermion singlets Aq, iV^, N2h (each of them having two degrees of freedom), with respective 
masses Mi, M2 — fJ-2, M2 + ^2 and Yukawa couplings given by the Lagrangian 



C Y = — A Q i h) P R N^ a - \ a2 P R N2h ^ N2l £ a + h.c. . (3.1) 

Recall that the parameters A a i violate lepton number and hence A a i <C X a 2- In turn this 
implies that the CP asymmetry in Aq decays is the dominant one. As a first approximation, 
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in the Eq. (3.1) we have neglected the L- violating couplings X' a2 <C \ a 2, because we expect 
that their contributions to the CP asymmetries and washouts are negligible. We have 
checked indeed that this is the case in the parameter space region relevant for leptogenesis 
that can also accommodate the observed light neutrino masses (see Sec. 4). As we will 
explain below, it is convenient to take Mi < M2 in order to obtain the lowest energy scale 
for leptogenesis within this framework. 

The amount of leptons and antileptons present in the Universe can be described by 
density matrices in flavour space. The evolution equations of these density matrices take 
the simplest form in the basis that diagonalize them, which are determined by the fastest 
interactions in flavour space (if there is a hierarchy among the different interactions). 
Let us suppose that the fastest interactions during Ai-leptogenesis are the A^-Yukawa 
interactions 1 . In this case the BE would be diagonal in a basis (£2, £2±, ^2_i_), with £ 2 ± 
and £*2± two orthogonal lepton flavour states perpendicular to the A^-decay eigenstate £2- 
But e(Ai -)• £ 2 ±h,£'2±h) = (with e(Ai -)• £ a h) the CP asymmetry in the decay of Aq 
into £ a h), and since the total CP asymmetry is null, then also e(Aq — > £2^1) = 0, with the 
result that no lepton asymmetry would be produced. Moreover, no asymmetry is generated 
when the Ai-Yukawa interactions are the fastest ones, because in this case the BE would be 
diagonal in the basis (£i,£i±,£'i±), and since the CP asymmetry in the Ai-decay eigenstate 
e(Ni — > £\K) coincides with the total CP asymmetry, it also vanishes. So we demand that 
the couplings of Aq and A2 be small enough, such that the Yukawa interactions of the 
r are the dominant ones. This implies that the BE are diagonal in the orthogonal basis 
(£ T , £ T ±, £' T ±), with £ T ± and £' T ± being determined by the fastest interaction acting in the 
plane perpendicular to £ T . For simplicity we take the Aq and A2 decay eigenstates to be 
perpendicular to £ e , in which case the aforementioned basis is £ T ,£^,£ e . Since no asymmetry 
is generated in the £ e flavour, two BE will be enough to determine the evolution of the 
lepton asymmetry and consequently the flavour indices a and (3 run over the species n and 
r. Later we will comment on the more general case with three non-null flavour asymmetries. 

Besides the charged lepton Yukawa interactions, the most relevant processes are the 
decays and inverse decays of the different heavy neutrinos, and the L-conserving but La- 
violating scatterings tph — > £ a h, £ph — > £ a h, and hh — > £ a £/3, hereafter called generically 
flavour changing interactions (FCI). We do not consider finite temperature corrections 
to the particle masses and couplings [18], moreover we also neglect spectator processes 
during leptogenesis and the asymmetry developed among the degrees of freedom of the 
Higgs [47, 48], as well as AL = 1 scatterings [49, 50] 2 . The relevant set of BE for the case 



1 When the effects of N 2 t and N 2 h simply add up we will refer to these states generically as N 2 . In 
the models considered here the only situation which requires a differentiated treatment of the degrees of 
freedom associated to N 2 is when fj, 2 <Si ^n 21 2h , because a density asymmetry can develop between the 
states produced by Ih and £h (see below). 

2 We have checked that the inclusion of spectator processes and AL = 1 scatterings modifies the results 
by at most a few tens of percent. 
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> rv 2li2h is 





= _1 1 


<Y Nl 


dz 


sHz \ 




dY Aa 






dz 


sHz 









1 7ft, (3-2) 



> , (3-3) 



where z = M±/T (with T the temperature), is the Hubble expansion rate, Yx = n^/s 
is the number density of a single degree of freedom of the particle species X normalized to 
the entropy density s, yx = (Yx — Y x )/Y x q is the asymmetry density normalized to the 
equilibrium density Y x 9 , and Y& a =Yb/3 — Yi a , with Yb the baryon asymmetry and Yz, a 
the lepton asymmetry in the flavour a. Since we are neglecting spectator processes during 
leptog enesis, the asymmetry in the lepton doublets £ a can be expressed in terms of Y& a 
using the simple relation Y<\ a = —Yi a = —2y( a Y^ [48]. The sphaleron processes convert 
part of the final B — L asymmetry (Yg_ L = y£ ) into a baryon asymmetry (Yg) and 
the relation between these quantities is taken to be Yg = 28/79 Y^_ L [51]. In Eq. (3.3) we 
have also introduced the notation 7°'j'"' = 7(0, b, . . . — > c, d, . . . ) for the reaction density of 
the process a, 6, ...—>■ c, c/, ... , the prime in 7^ indicating that the on-shell contribution 
has to be subtracted, and 7^ = J2 a 7i3i+7^V Finally notice that the lepton asymmetries 
generated by N2 decays, e a 2, have not been taken into account, given that e a 2 <C e a \. 

Instead, if ^2 "C rjv 2i2h then N21 and iV^h combine to form a Dirac neutrino N2 = 
(N2h + iN2i)/V2, and therefore there is an asymmetry generated among the degrees of 
freedom of N2 which has to be taken into account [46]. An appropriate set of BE for this 
case is 

~d7~ = 7Wz yy^ ~ J 1 ' (3 - 4) 

%^~E7S^ 2 -|/d, (3-5) 



dz 

dY Aa -1 J / Y Nl 
dz sHz 1 \ Y e N \ 



/3^a J 

The reaction densities for the FCI are obtained integrating the corresponding cross sec- 

T r°° ( fs\ 

tion ctfcI) 7fci = 7 / ds s^ctfci-^i ( — ) > with the reduced cross sections cxfci = 

64vr 4 7 V T J 

2so"fci- In order to give general expressions for these cross sections we take the FCI to be 
mediated by any number of Dirac or Majorana neutrinos TVj (i = 1, 2, . . . ) with masses Mj, 
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whose Yukawa couplings with the lepton doublets £ a are called h a i. We find 3 : 
a(tph ->■ £ a h) 



1 ^-A^-iJWirO^-JWj + tM^-) 

^ 2^ h ai n l3i tl ajri l3j 



1,3 



Mf) + M 2 rf ( S -M|) z + M|r2 



i,j 1 3 V 

i+3 

1 ST^u* u u u* fi S + M i 
+ 2^2^ h <xiHi h <*3 h Vj I In — ^ 



2i s + M? 



ill 2 



M 2 In ■ 



+ M 2 



a(hh — > t a i 



— ^ h ai h pi h aj h^ -1 + , m2 _ M?) ^ - 



ill 2 



M 2 



M 2 (, + M 2 ) in , + M 2 
s(M 2 - M 2 ) M 2 

a + 2Mf s + Ml 
s n Mf 



Moreover, the subtracted reaction density is given by 

tph' _ iph sr^ tfih-, 



s + Mf 

~my 



+ 



(3.7) 



(3. 



with Br(iVj — > £ a h) the branching ratio of the process iVj — > £ a h. 

For the models we are considering the dominant contribution to the FCI comes from 
the exchange of a neutrino with mass M2 and Yukawa couplings h a 2 = X a 2 4 - 



4. Results 

We will find the lower bound for M\ under the requirement that leptogenesis be successful in 
these scenarios. The relevant parameters for leptogenesis are Mi, M2/M1, (A' A)n, (A^A)22, 
the projectors K a i = A a jA*j/(A^A)ji (a = n,r;i = 1,2), and \i2- Next we comment on the 
role of these parameters. 

• M\: For a fix value of the Yukawa couplings (and in particular of the CP asymmetry 
when M 2 /Mi is kept fixed) the lower is M 1; the stronger the washouts become (be- 
cause the expansion rate becomes slower with decreasing T) . This is why M\ cannot 
be very low. 

3 Our expressions for these cross sections differ from the ones used in [39], which were taken from [26]. 

4 Note that the contribution to the FCI of the two Majorana neutrinos N 2 i , N 2 h with the Yukawa couplings 
of Eq. (3.1) is the same as the contribution of one Dirac neutrino with mass M2 and Yukawa couplings X a2 , 
in the limit fj,2 <C M2. 
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• M 2 /Mi : The hierarchy among the masses appears in the CP asymmetry, with e Q i oc 
(M 2 /Mi)~ 2 for M 2 S> Mi, and in the FCI when these are mediated mainly by iV 2 
(as always happens in the interesting cases for this work), jfci(T) oc (M 2 /Mi) -4 for 
T ~ Mi < M 2 . Also note that if M 2 /Mi < 20 the inverse decays of iV 2 could erase 
part of the asymmetry generated during Ni leptogenesis and therefore they should 
be included in the BE. 

• (A^A)n: The intensity of the washouts due to processes involving Ni is determined by 
the effective mass rh\ = (A^A)nw 2 /Mi. In order to make full use of the -mandatory- 
flavour effects we take fhi > m* ~ 10~ 3 eV 5 . However, since (A^A)n violates lepton 
number it cannot be very large in the scenarios we are considering. 

• (A^A) 22 : The CP asymmetry e a i is directly proportional to (A^A) 22 , hence this pa- 
rameter should be taken as large as possible. There are two reasons that prevent 
it from being very large. One is that the washouts due to processes involving iV 2 
increase with (A^A) 22 . The other is that, as explained at the beginning of Sec. 3, the 
Yukawa interactions of A 2 must be slower than those of the r. Below we explain in 
more detail how this constraint has been dealt with. 

• K a i\ The flavoured CP asymmetries depend on the square root of the projectors 
while the washouts in the flavour "a" depend linearly on K a i, hence the washouts 
decrease faster than the CP asymmetries with decreasing projectors. This fact must 
be taken into account in order to maximize the production of lepton asymmetry. Note 
that since J2 a ^ai = 1 and we are taking K e i = 0, there are only two independent 
projectors, which can be chosen as K^i and K^- 

• [12 : Since we are not considering the case /U 2 ~ T n 21 2h , /i 2 only enters as a discrete pa- 
rameter, the baryon asymmetry taking one /x 2 -independent value when jit 2 3> Tjy 2l 2h 
and another when ^ 2 <C Fjy 2l 2h . 

Given that the hierarchy M 2 /Mi is an interesting and crucial parameter, we have 
determined the minimum value of Mi compatible with successful leptogenesis as a function 
of Mi/Mi, maximizing the final baryon asymmetry over the remaining parameters, i.e. 
over the relevant combinations of Yukawa couplings: (ATA)n, (A' A) 22 , Kn\, and K^- To 
obtain the baryon asymmetry we have solved numerically the appropriate set of BE, and 
to get successful leptogenesis we have required Yb = 8.75 x 1CT 11 [52]. The result is 
represented with the thick continuous curves in Fig. 1, the red line corresponding to the 
case H2 > ^N 2 i.2h an d the green one to /i 2 <C ^N 2l 2h . 

As can be seen it is possible to have neutrino masses as low as Mi ~ 10 6 GeV, i.e. 
around three orders of magnitude below the lower bound for the standard case of type I 
seesaw with hierarchical heavy neutrinos. An interesting consequence is that this value 

5 The quantity m* is the equilibrium mass, which is defined by the condition H ^} Ml ^ = so that 

m * ~ ttf 71 " 5 ^ 2 \/9*sm-t— ~ 1,08 x 10~ 3 eV (q*sm is the number of SM relativistic degrees of freedom at 
temperature T and m p i is the Planck mass). 
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M2/M1 

Figure 1: Lowest value of Mi yielding successful leptogenesis as a function of M2/M1. The red 
curves are for the case /12 3> Tn 2 i 2h and the green ones for fjQ "C IV 2( 2h - The thick continuous 
curves give the physically correct bound, while the thin dashed ones show the result that would be 
obtained if the Yukawa couplings of N2 were allowed to take values as large as 1 for all values of 
M 2 /Mx. 



can be compatible with the upper bound on the reheating temperature required to avoid 
the gravitino problem in SUGRA models. Moreover, M\ values around 10 6 GeV can be 
achieved for a wide range of N2 masses and also for different values of the Yukawa couplings. 
In particular, the maximum value of the baryon asymmetry is obtained for rh\ ~ 10~ 2 eV 
and Knx (or K T \) ~ 0.1, but it does not change much if rh\ is varied within the range 
5 x 10~ 3 eV < fh\ < 1CT 1 eV as long as the projectors are adjusted in order to have 
weak washouts in one flavour and strong washouts in the other, e.g. K^rhx < and 
K T iihi > (5 — 10)771*. Regarding the Yukawa couplings of N2, (A^A)22 lies approximately 
between 6 0.01 and 1, while for intermediate to large M2 masses, M2 > 5Mi, the smallest 
projector can take values as large as 0.1 to 0.5, without changing the bound shown in 
Fig. 1 by more than a factor 2. Instead, for M2 < 5Mi, some hierarchy among the Yukawa 
couplings of A^2 is required, namely A M 2/A r 2 (or A r 2/A M 2) ~ (1 — 3) x 10~ 2 to achieve the 

6 In the analysis presented here we have restricted | Aq. 2 1 < 1, but we have checked that the results do not 
change significantly when allowing these Yukawa couplings to take somewhat larger values (but below the 
perturbative bound). 
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values shown in Fig. 1. 

An important issue for obtaining the bound on Mi is to determine how large the 
Yukawa couplings of A^2 can be without violating the condition that the rates of processes 
involving N 2 be slower than the rates of the r- Yukawa interactions 7 . For reference purposes 
we plot in Fig. 2 the ratio between the rates of these interactions, distinguishing between 
processes with a real or virtual N 2 , for M 2 = W 7 GeV and v / (AtA) 2 2 = 0.01 ~ h T (the 
rates of the r- Yukawa interactions have been taken from [53, 54]). On one hand, note that 
the r- Yukawa interactions tend to become dominant over the N 2 ones as the temperature 
decreases. On the other hand, for each set of parameters there exists a temperature Tj n 
such that the final baryon asymmetry (Yg) does not depend on what happens at T > Ti n 
(because there are strong washouts for all the interesting regions of the parameter space). 
Then the largest possible value for (A^A) 2 2 has been set by requiring that the r- Yukawa 
interactions be faster than the N 2 ones for all temperatures below T{ n . In practice we have 
determined T{ n finding the lowest value of T such that Yg does not change by more than 
10% when the initial conditions at T, namely Yjv^T) and Y& a (T), are varied. The value 
of T in as a function of M 2 /Mi is shown in Fig. 3, reinterpreted as an approximate lower 
bound for the reheating temperature. Also note that the independence from the initial 
conditions is an interesting feature by itself, since the models become more predictive. 

For comparison we have also plotted in Fig. 1 the -wrong- bound that would be obtained 
if (A^A) 2 2 were allowed to be as large as 1, ignoring the above discussion. It is clear that as 
M 2 approaches Mi the requirement of an upper bound for (A^A)22 becomes very relevant. 
This requirement is also the reason why it is convenient to take M 2 > Mi, because in the 
opposite case (when the lepton asymmetry is produced by the decay of the next-to-lightest 
singlet neutrino and that is usually called U N 2 leptogenesis" ) the bound on (A^A)22 would 
not allow for large values of the CP asymmetry. 

Another feature apparent in Fig. 1 is the constant behaviour of the bound on Mi for 
large values of M 2 jM\. This can be understood quite easily as follows. The asymmetry 
Yg is proportional to the CP asymmetry in N\ decays, e M i = —e T \, which in turn satisfies 
e M i oc (Mi/M2) 2 (A At 2A T 2). For M2 ~> Mi, Yg increases with (A^A T 2) up to a certain value 
(A jU 2A T 2)max for which the FCI become important, equilibrating the asymmetries generated 
in the two flavours with the consequent decrease of the baryon asymmetry [40, 41]. At 
the temperatures relevant for leptogenesis, T ~ Mi <C M2, the rates of the FCI scale 
as 7fci oc (Mi/M2) 4 (A At 2A T 2) 2 - Hence we see that if (M1/M2) is decreased by a factor a, 
(A A t2A T 2)max can be increased by a factor a 2 keeping constant the rates of the FCI. Therefore 
(A M 2A T 2)max oc {Mi/M 2 )~ 2 and hence the maximum value of Yg, being proportional to 
(M1/M2) 2 (A^2A T 2)max, becomes independent of M 2 jM\. 

As mentioned in Sec. 3, for simplicity we have taken l e to be perpendicular to N\ and 
N 2 , so that only two flavour asymmetries are generated. We have checked that in the more 
general three flavour case it is possible to lower the bound on Mi by a factor up to almost 
4 with respect to the two flavour case. This is due to the combination of two effects. On 
one hand, Yg can be almost twice as large compared to the two flavour case (this happens 

7 There are not relevant experimental bounds on the Yukawa couplings of Ni for the masses of heavy 
neutrinos we are considering. 
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M 2 = 10 7 GeV, [V r k) 22 = W A 




10000 100000 1e+06 1e+07 1e+08 1e+09 
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Figure 2: Comparison between the rates of the Yukawa interactions of N 2 and r as a function of 
the temperature. The red line gives Jt/jn 2 and the green one j T / |tsfci|, where j T is the rate of 
the Yukawa interactions of the r, jn 2 is the sum of the rates of processes involving a real iV 2 (here 
we have included the decay and scatterings with the top quark summing over all lepton flavours), 
and 7efci is basically the sum of the rates of the FCI mediated by an off-shell N 2 , 7efci = 
(jth ' + 7t T Jh + Ihh ) /{K&Kri) . We have taken M 2 = 10 7 GeV and (AtA) 22 = 10" 4 ~ h* (with 
h T the Yukawa coupling of the r). Note that the interactions involving a real N 2 scale as (X'X)22> 
while the FCI are proportional to (A'A) 2 2- The corresponding curves for other values of M 2 can 
be obtained simply by making the appropriate translation along the T-axis. The spiky shape of 
the green curve for 0.1M 2 < T < M 2 is due to the substraction of the on-shell contribution to the 
process i T h — > i^h. 

when K a ifhi < for two different flavours). On the other hand, it can be shown that 
for a given value of Mi/Mi, the maximum value of Yg is proportional to ~ \/M\ in the 
relevant region of the parameter space. Hence an increase in Yg by a factor 2 leads to a 
decrease in the lower bound of M\ by a factor 4. 

Finally let us comment on the relation between the parameters defined above and the 
light neutrino masses. To lowest order in the L- violating parameters, the light neutrino 
masses m u are given by 

(m„W ~ A al ^-A^ + (X' a2 - ^A Q2 )^-A^ 2 + A a2 ^-(A^ 2 - ^\ p2 ) (4.1) 
where v = (h) = 174 GeV is the vev of the Higgs field. 
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Figure 3: The lower bound on the reheating temperature as a function of MijM\. The red curve 
is for the case /12 3> rV 2i 2h and the green one for [12 <C IV 2i 2h . 



We have seen that M\ is minimized for values of rh\ = (A^ X)nv 2 jM\ in the range 
1CP 2 eV < rh\ < 10 _1 eV, therefore the contribution of N\ to light neutrino masses is 
expected to be of the same order, barring cancellations due to phases. Imposing that such 
contributions are of order m a tm ~ 0.05 eV, we get A a i ~ 10" 5 - lO -4 . 

To reproduce the observed solar and atmospheric mass scales, at least one of the 
contributions from N2 in (4.1) should be also of order 10 -2 eV; for the parameters that 
minimize Mi, this requirement leads to ~ 10~ 8 — 10~ 6 , independently of the ratio 

M 2 /Mi. Moreover, T N2l2h /M 2 ~ 5 x (10^ 4 - 10" 2 ), therefore typically fj, 2 <C ? N2l2h . 
However, for M\ > 5 x 10 6 GeV, and still not considering large fine tunings related to 
phase cancellations, smaller values of \ a 2 can lead to successful leptogenesis, and in this 
region it is possible to have 112 ^ ^n 2 i 2h- 

With respect to the L- violating parameters A^, 2) their contribution to the masses of 
the light neutrinos is m v ~ rn a t m typically for X' a2 ~ 10~ 8 — 1CT 7 . We have checked for 
consistency that these small values give a negligible contribution to leptogenesis (more 
specifically to the CP asymmetries and washouts). 
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5. Conclusions 



We have studied leptogenesis in the framework of the seesaw mechanism with small viola- 
tion of B—L. If B—L is only slightly broken, then either the heavy neutrinos which generate 
the baryon asymmetry are almost degenerate and combine to form a quasi-Dirac fermion, 
which can have large lepton number conserving Yukawa couplings, or it is a Major ana 
fermion with small lepton number violating Yukawa couplings. 

In both cases there are interesting consequences for leptogenesis: in the first one, 
the strong degeneracy of the heavy neutrinos leads to a resonant enhancement of the CP 
asymmetry, avoiding the DI bound on Mi which applies to hierarchical SM singlets. In 
the second case, the L-conserving part of the flavoured CP asymmetries in N\ decays can 
be much larger than the L-violating one, since the former is not linked to light neutrino 
masses and also escapes the DI bound, even if the heavy neutrinos are hierarchical. As a 
consequence, in models with almost conserved B — L successful leptogenesis may be possible 
at lower temperatures than in the standard seesaw, alleviating the gravitino problem in 
supersymmetric scenarios. 

In this paper we have focused on the second possibility, i.e., we have not resorted to the 
resonant enhancement of the CP asymmetry. Since the L-conserving contributions to the 
flavoured CP asymmetries cancel in the total CP asymmetry, e» = ^2 a ej a , it is mandatory 
that flavour effects are at work for these terms to have an impact in leptogenesis. 

We have exhaustively scanned the parameter space of seesaw models with almost con- 
served B — L, in which the L-conserving piece of the flavoured CP asymmetries dominates 
(purely flavoured leptogenesis). We have found that the largest baryon asymmetry is gen- 
erated by the lightest SM singlet, and it rapidly decreases if |M2 — M±\ < 2M\ (but far 
from the resonance region, \M% — M\\ ~ T2), being N2 the next-to-lightest heavy neutrino; 
thus we have restricted the heavy neutrino masses to the region M2 > 2M±. The relevant 
parameters are Mi, M2/M1, (A^A)n, (A^A)22 and the flavour projectors K a i. The same 
A^-Yukawa couplings which enhance the L-conserving CP asymmetries induce large FCI 
mediated by 7V 2 , which tend to equilibrate the asymmetries in the different lepton flavours 
diminishing the total lepton asymmetry, especially if both N\ and N2 have similar masses 
and are simultaneously present in the thermal bath. Therefore, for each value of M2/M1 we 
have determined the minimum Mi compatible with successful leptogenesis, maximizing the 
final baryon asymmetry over the remaining parameters. We have solved numerically the 
relevant set of BE, including decays and inverse decays of the two singlet neutrino species, 
as well as the FCI. We have considered two possibilities: N2 is a pseudo-Dirac fermion (i.e., 
two Majorana neutrinos with masses M2 =t fi2, 1^2 *C M2) and N2 is approximately Dirac, 
in which case an asymmetry between N2 and N2 is generated and should be taken into 
account. In both cases we have found that leptogenesis is possible for Mi > 10 6 GeV, as 
long as M2/M1 > 5 and (A^A)22 ~ 0.01 — 1 (see Fig. 1). So purely flavoured leptogenesis in 
seesaw models with slightly broken B — L provides a solution to the conflict between the 
upper bound on Trjj required to solve the gravitino problem of supersymmetric scenarios 
and the lower bound on Trh needed for successful thermal leptogenesis. However, such 
heavy neutrinos are far outside the reach of present and near future colliders and do not 
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lead to observable lepton flavour violation in non-supersymmetric frameworks. 
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